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Radiative Transitions in Charmonium from Lattice QCD 
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Radiative transitions between charmonium states offer an insight into the internal structure of 
heavy-quark bound states within QCD. We compute, for the first time within lattice QCD, the 
transition form-factors of various multipolarities between the lightest few charmonium states. In 
addition, we compute the experimentally unobservable, but physically interesting vector form-factors 
of the rjc, J/tp and x^o- 

To this end we apply an ambitious combination of lattice techniques, computing three-point func- 
tions with heavy domain wall fermions on an anisotropic lattice within the quenched approximation. 
With an anisotropy £ = 3 at a s ~ 0.1 fm we find a reasonable gross spectrum and a hyperfine split- 
ting ~ 90MeV, which compares favourably with other improved actions. 

In general, after extrapolation of lattice data at non-zero Q 2 to the photopoint, our results agree 
within errors with all well measured experimental values. Furthermore, results are compared with 
the expectations of simple quark models where we find that many features are in agreement; beyond 
this we propose the possibility of constraining such models using our extracted values of physically 
unobservable quantities such as the J/tp quadrupole moment. 

We conclude that our methods are successful and propose to apply them to the problem of 
radiative transitions involving hybrid mesons, with the eventual goal of predicting hybrid meson 
photoproduction rates at the GlueX experiment. 

PACS numbers: 12.38.Gc, 12.39.Pn, 12.39.Jh, 12.40.Yx, 13.20.Gd, 13.40.Gp, 14.40.Gx 



I. INTRODUCTION 

Charmonium occupies a valuable intermediate position 
within QCD, being neither in the purely non-rclativistic 
regime nor the regime where chiral symmetry breaking 
dominates the physics. This makes it a relatively clean 
system in which to study non-perturbative QCD dynam- 
ics, and indeed QCD-inspired quark-potential models as 
well as lattice QCD have been rather successful in de- 
scribing the observed features of the spectrum^]. How- 
ever, charmonium cannot be considered to be completely 
understood; as an example, in recent years a number of 
new charmonium resonances have been claimed in exper- 
iment (see @ for a review), several of which cannot be 
easily reconciled with the predictions of simple quark- 
potential models. 

Unlike the light quark sector, in charmonium the light- 
est state for most J 's lies below the threshold for OZI- 
allowed decay and consequently these states are rather 
narrow. These states have been the subject of many cal- 
culations in lattice QCD which generally reproduce the 
gross structure of the spectrum, but are unable to ac- 
count for all of the detailed fine structure (such as the 
117 MeV J/tp — i] c splitting), owing to some combina- 
tion of the approximations inherent in the computations, 
which can include the finite lattice spacing, quenching 
and lack of disconnected diagrams. 

Masses are not the only well measured charmonium 
observables. Because of the small total width of these 
lightest few states, radiative transitions between them 
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constitute large branching fractions and have been mea- 
sured experimentally by a number of groups 0-Q These 
quantities have been studied within quark-potential mod- 
els (and latterly EFT approaches like pNRQCD) where 
they are related to the overlap of meson wavefunctions 
with the photon current operator and as such are an in- 
sight into the internal structure of these states. QCD 
sum rules have also been applied with some success |5( . 
No study of radiative transitions in charmonium has yet 
been performed using lattice QCD - it is this situation 
that we remedy in this paper. 

This study in the charmonium sector is an ideal test- 
bed for our eventual aim of computing photocouplings 
in the light-quark sector, in particular the coupling be- 
tween a conventional meson, a photon and a hybrid me- 
son. Such a coupling drives the photoproduction mech- 
anism proposed by the GlueX collaboration for their ex- 
periment in Hall D of the upgraded CEBAF at Jefferson 
Lab@. Flux-tube model calculations 0, @ suggest that 
these couplings are not small and that we may expect co- 
pious production of hybrid mesons, but the assumptions 
underlying such a model need to be tested in a frame- 
work closer to QCD. This paper will focus on transitions 
between conventional charmonium mesons, and we will 
proceed to transitions to hybrid mesons in a subsequent 
work. 

In section[n]we outline the computational details of our 
lattice calculation and in section ITTT1 we present the char- 
monium spectrum so obtained. In section ll*Vl we explain 
how three-point functions extracted from the computa- 
tion can be converted into multipole form-factors and in 
section^we present lattice estimates for the vector form- 
factors of the lightest three charmonium states and com- 
pare with what one would expect on the basis of simple 



2 



quark models. In section lvTl we consider radiative transi- 
tions and discuss the extrapolation from the non-zero Q 2 
accessible on our finite- volume lattice to the relevant case 
of an on-shell photon. Finally in section lvTTl wc conclude, 
relegating the details of Lorentz covariant multipole de- 
compositions and scale setting on anisotropic lattices to 
appendices. 

II. COMPUTATIONAL DETAILS 

The computations were performed in the quenched ap- 
proximation to QCD, using the Chroma software sys- 
tem 0. We employed 300 configurations on a 12 3 x 48 
lattice generated using an anisotropic Wilson gauge ac- 
tion |lOj . with a renormalized anisotropy £ = a s /a t = 3. 
The temporal lattice spacing obtained from the static 
quark-antiquark potential is a^ 1 = 6.05(1) GcV. 

Anisotropic lattices as applied to charmonium exploit 
the fact that while the quark mass scale demands a 
cut-off above ~ 1.5 GeV, the internal three-momentum 
scale is typically much lower, ~ 500 McV. On our lat- 
tice, we can have both m c at and \p\a s reasonably small 
and a spatial length > 1 fm without requiring very 
many spatial lattice sites. We work on only one vol- 
ume, L s sa 1.2 fm; previous charmonium spectrum stud- 
ies indicate that there are no significant finite volume 
effects for lattices of this size or larger [TlllT2l |. With this 
volume and anisotropy the three-momenta accessible are 

The quark propagators were computed using an 
anisotropic version of the domain-wall fermion (DWF) 
action [l3| with a domain-wall height ajM = 1.7, a 
fifth dimensional extent L5 = 16, and a quark mass 
a t m q . The conventions used for the action are defined in 

Ref. 0| . In terms of dimensionless variables = a^ 2 ip 
and — a^Wp, where ao = a t and = a s , k = 1, 2, 3. 
The kernel of the domain action is the anisotropic Wilson 
fermion action 0, : 

a t Q = -a t M + v t W tli + ^ V W h lk, (1) 

where the Wilson operator 

W» = V M - ^7mA m • 

The parameter £0 i n Eq.^is the bare anisotropy which is 
determined in the gluonic sector so as to yield the desired 
renormalized anisotropy £. The remaining parameters of 
the action, m q ,v s and v t represent the quark mass, and 
the renormalization of the couplings in the spatial and 
temporal directions, respectively. The parameters v s and 
v t are not independent; we will fix v t = 1 and then tune 
(m q , v s ) so as to yield the desired quark mass and such 
that the speed of light obtained from the meson disper- 
sion relations be 1.0, as discussed in the next section. 

While this computation is performed at only a single 
value of the lattice spacing and hence no attempt can be 



made to determine the lattice spacing dependence of the 
results, an important benefit to the use of the domain- 
wall fermion action is that it is automatically improved 
to C(a 2 ). 

We anticipate that the quenched approximation will 
not be a serious impediment to this study. In the sub- 
DZ)-threshold charmonium system we expect (on the ba- 
sis of previous lattice studies) that the dominant effect of 
light quarks will be to modify the running of the strong 
coupling. In the quenched approximation the QCD /3 
function is not properly reproduced, thus if we set our 
lattice scale using a long distance dominated quantity 
(such as the IP — 15 splitting or some intermediate dis- 
tance in the static quark-antiquark potential) then we 
will have a weakened coupling at short distances com- 
pared to QCD. Hence within the quenched approxima- 
tion we cannot expect to correctly describe such short- 
distancc-dominatcd quantities as the hypcrfinc splitting, 
or meson decay constants |l7l ITsI ]. Our principal inter- 
est is with radiative transitions, which are dominated by 
long-distance wavefunction overlaps, and hence should 
not be considerably distorted by the quenched approxi- 
mation. 



III. SPECTRUM 

Charmonium masses and interpolating field overlap 
factors were extracted from fits to the connected two- 
point functions 1 , T^(p;t) = 

Z (N) Z (N)* 

£ e-^rw, t) #>(o, o)) = J2 1 2E( n) e ' E<N 

x N 

(2) 

The fact we used local fermion bilincars as interpolating 
fields limits us to the J PC listed in Table [I] 

In the previous section we discussed the need to have 
a rather fine lattice spacing to ensure that charm quarks 
are not "cut-off" ; a related problem arises from the two- 
scaledncss of the charmonium system. The scale which 
dictates our choice of lattice spacing is the typical mass of 
a charmonium state ~ 3GeV, while the level spacing e.g. 
J/tp — i// ~ 600MeV, is much smaller. This second scale 
indicates how fast, relatively, the excited state exponen- 
tials in eqn J5J) fall off with t, which we see will be rather 
slow, so that it is possible at finite times to get contribu- 
tions from many excited states. This is not the case with 
light quarks where the level splitting is of the same size 
as the ground state mass (e.g. p(770), p(1450)). Because 
many excited states can be contributing during the finite 
temporal-extent of our lattice, we are unlikely to reach a 
plateau and as such extracting the ground-state mass is 
non-trivial. 



1 Disconnected diagrams will be discussed in section II VI 
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We can reduce this problem somewhat if we improve 
the overlap of the interpolating field on to the ground 
state by smearing it over space. We do so using a gauge 
invariant cubic approximation to a rotationally symmet- 
ric Gaussian, 
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With the right smearing radius this will resemble the 
ground state wavefunction and maximise the overlap, 
while the excited states, which have radial nodes, will 
have a decreased overlap. 

We computed two-point functions with smearing radii 
(7 = 0.0,3.6,4.0. For the case a = 3.6, we attempted 
to reduce gauge fluctuations by smearing the gauge links 
entering into cqn. [31 using the gauge- invariant procedure 
(APE) introduced in rcf. For a = 3.6,4.0 we used 

TV = 32, 50 and 5, APE smearing iterations. 

We display effective mass plots in figure ^ where the 
oscillations in the first few timeslices are an artifact of the 
non-locality of the DWF action in four dimensions. The 
(7 = 3.6 smeared-smeared data has reasonable plateaux, 
at least for the rj c , if) and Xco- 

Our fitting procedure is as follows. Given a set of in- 
terpolating operators {O;}, we can construct the matrix 
of time-sliced correlators 

rf(p;t) = ]Te-^(O^,i)Oj(0)) 



= E 



AN), 



N 



2E(N) 



(4) 



where 



zf\p) = (Q\O l \N{p)). 



(2) 

If we had computed every element of the , then we 
could apply a variational method to diagonalise the sys- 
tem. However, in general we have an incomplete ma- 
trix, with several sink operators for each source oper- 
ator. Thus we employ a "factorizing" fit to this incom- 
plete matrix of correlators, in which the fitted parameters 
1, . . . , N op , i = 1, . . . , N cxp }, where 



are {E®,zf ] : j 



there are N op operators used in the fit, and we include 
iVcxp states. We assume that the residues of the fit fac- 
torize according to eqn. 0] and that as demonstrated in 
Appendix El the Z's are real. In practice we include the 
three smearing combinations shown in figure ^ for each 
J PC and where possible fit to the ground state plus one 
excited state. For the xa and h c the noise on the data 
allowed for only inclusion of the ground state in the fit. 

In order to propagate statistical error due to the finite 
number of gauge field configurations through our calcula- 
tions we adopt a single elimination jackknifc procedure. 



N 
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FIG. 1: Effective mass plots (at p = (000) using local inter- 
polating fields. Dashed lines and rightmost data point show 
the ground-state mass obtained from the multi-correlator fit. 
(a) smeared(<r = 4.0)-local (b) smeared(<r = 3.6)-local (c) 
smeared(o" = 3.6)-smeared(cr = 3.6) 



For the two-point function fits we remove from the en- 
semble one configuration and average the remaining two- 
point functions. The factorising fit is performed on this 
average yielding a value for each of the fit parameters 
(Z,E). This is repeated for each configuration yielding 
an ensemble for the fit parameters over configurations. 
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These are saved for later use in the three-point function 
calculations. 

Fits are shown in Figure [5] and Table [I] 
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We extract a J/ip — i] c hypcrfmc splitting of 97(6)McV. 
This is probably 10% larger than it would be if we had 
correctly tuned the quark mass, assuming the hyperfine 
splitting scales like m~ 2 ; with this correction we are in 
rather good agreement with the result 88(4)MeV from 
the closely related overlap fermion action on an isotropic 
lattice at Wilson ft = 6. 3^^. These two results are sig- 
nificantly larger than typical results using other improved 
actions such as clover. Assuming that there is very little 
change in scaling to the continuum this leaves a deviation 
from experiment of ~ 30MeV, which is consistent with 
the sort of deviations suggested by un a ue nching [l7l Elj 
and inclusion of disconnected diagrams [221 |23| . 

The excited state masses we are able to extract are 
clearly too large in comparison with experiment. The 
same effect was observed in the anisotropic clover study 
of [24( when they set their scale using ro, and to a lesser 
extent in the studies of ^(| and [25| . One possible expla- 
nation may be that, as discussed above, several excited 
states are contributing at intermediate times and that 
our fitting is not able to fully resolve them such that we 
end up with one effective excited state "mocking up" the 
effect of several 2 . Including a second excited state into 
the fit hypothesis did not produce stable fits. 

The multi-correlator fits were performed on the data 
with many p to determine the dispersion relation for 
our meson states. This is particularly relevant on an 
anisotropic lattice where parameters in the fermion ac- 
tion must be tuned to match the gauge sector anisotropy 
and set the speed of light c 2 = E ^l^" 1 equal to one. 
The results are shown in figure |3| where we observe that 
our tuning was good to the 3% level onc= 1.032(9). 



2: Data(cr = 3.6 — a = 3.6)-fit percentage deviation, 
= ioo data-fit _ Black points indicate the contribution of 

d at a, 1 



FIG. 

A (2) 

the ground-state to the fit, the red points include both the 
ground and first excited states, (a) 770(75), V — (000), (b) 
^(75), P = (100), (c) ^(Tfc), p= (000), (d) 7>( 7fe ), p = (100), 
(e) Xoo(l), P= (000), (f) Xco(l), P= (100), (g) Xci(7fc7s), 
p = (000), (h) M7fc74), P= (°°°) 
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TABLE I: Local interpolating fields, tpFtp and extracted char- 
monium spectrum, using the determination of the lattice 
spacing as described in the text. 




FIG. 3: Speed of light extracted from meson dispersion rela- 



tion: c 



The ground state masses extracted are observed to be 
systematically around 5% lighter than the experimental 
masses - this is a result of imperfect tuning of the quark 
mass parameter; for this first study of radiative transi- 
tions we did not attempt a more accurate tuning. 



2 Consider in particular the vector channel where the second ex- 
cited state i/i(3770) is very near to the first ?/>(3686), although its 
3 £>l nature may lead to a reduced overlap with the local inter- 
polating field 
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A. Charmonium decay constants 

Apart from the hyperfine splitting, the only other 
quantities we consider that might be significantly dis- 
torted by the quenched approximation are decay con- 
stants, which are related to the meson wavcfunction 
at the origin and are hence sensitive to short-distance 
physics. We discussed above that since we set the scale 
using long-distance dominated quantities, we expect dis- 
tortion of short-distance quantities owing to the incor- 
rect running of the QCD coupling. To be specific, the 
quenched coupling at short distances is too small and 
as such we expect the wavefunction at the origin to be 
depicted and hence the decay constants to be too small 
with respect to experiment. 

We define the rj c and ijj decay constants via 3 

(0\V»\i>((f)) = (0^(0)7^(0)1^(9^)) = mi> f^(q, A). 

(5) 

In the pseudoscalar (P = i/j"f 5 ip) case we can form the 
following object 



Zi°\e- m<0)t 



(48-t)) ' 



where the subscript a indicates that the interpolating 
field is smeared with radius a and where the superscript 
(0) indicates the ground state. The renormalisation of 
the local axial current, Zj±, will be discussed in the next 
section. Using (0), f Vc is equal to 



^ 0)+f vi 1} ~(o) 



-m (1) (48-t) 



-mWt 



-m(°>(48-i) 



(0) 



We fit the lattice data over the range t = 9 — 22 with this 
form using the zi 0,1 \ to' ' 1 ) found in the earlier two- 
point function fits and considering f^'^ to be free fit 
parameters - the result for the smearing choice a = 3.6 
is shown in figure 0Ja). 

We can perform an analogous fit to the vector two- 
point data, 



3^ 0) (e-" l(0)t + e- m(0) (48-*)) 



Urn +fj,m (0) (0)t 



( + e -m (1, (48-t) 



-m(°)(48-i) 



(7) 



the resulting fit with a = 3.6 being shown in figure Qfb). 
The renormalisation of the local vector current, Zy is 
discussed in the next section. 



Experimentally, the vector decay constants can be ex- 
tracted from the leptonic decay widths 4 , 



r(v 



f 2 

-) = ^2c?l± t 



47T 4 

Y9 ( 



so that, using the PDG averages^ we have fj/^ = 
411(7)MeV and f 4 ,> = 279(8)MeV. The ?y c decay con- 
stants are rather more difficult to measure; the only result 
available comes from B — > r\ c K [2^ and involves using a 
factorisation approximation to yield = 335(75)MeV. 
It is clear that within errors fj/^ ~ f ric as one would 
expect from the non-relativistic quark model in which 
the J/ip and the r\ c differ only by v/c suppressed spin- 
dependent terms. 




l — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — r 
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FIG. 4: Data points and fits (t = 9 — 22) - see equations @, 
Q. The black line indicates the ground state contribution, 
p°\ Upper plot - rjc, lower plot - ip. 



3 An alternative definition for the vector decay constant is often 
used in lattice simulations 



4 Note that we d o not include the explicit "QCD correction" factor 
(1 — -j^ois) l26l under the assumption that all gluon effects, even 
hard gluons close to the charm quark mass scale, are included in 
the path integral computed 



From our fits we extract (the scale setting on an 
anisotropic lattice is discussed in appendix Q 

\f Jfi ,\ = 399(4)MeV; |/„J = 429(4)(25)MeV 
|/^| = 143(81)MeV; \f v >J = 56(21)(3)McV 

The first error is statistical and the second, estimated 
systematic error, is of order 5% and accounts for the 
small violations of Zy = Za we find for our domain-wall 
fermions, see the following section for details. 

The ground state values are in reasonable agreement 
with the experimental extractions, perhaps indicating 
that even here, with a short-distance dominated quan- 
tity, the quenched approximation is only impacting at 
the few percent level. It is important however to recall 
that our charm quark mass is a little small, and it is 
possible that the discrepancy caused by this is approx- 
imately canceling the quenching discrepancy. A further 
calculation with a slightly larger quark mass would allow 
for an interpolation to the true charm quark mass and 
a more precise estimation; since our primary interest is 
with radiative transitions, this additional computation 
has not yet been carried out. 

The excited state decay constants extracted are some- 
what smaller than one might expect, however we recall 
that we were not very successful in extracting the exper- 
imental masses for these states from two-point functions. 
It is possible that the excited state exponential is absorb- 
ing the effect of several excited states - a more systematic 
study of interpolating fields with the intent of isolating a 
field with large overlap on to the first excited state would 
clarify the issue. 

IV. THREE-POINT FUNCTIONS 

The field theoretic quantity that most simply encodes 
radiative transitions is the three-point function: 

r /ri(p> »?;*/.*) 

= J2 e-^f s e +i ^(O f (x, t f ) t) O](0, 0)) (8) 

The Oij are interpolating fields chosen to have some 
overlap with whatever meson states we are interested in. 
The current insertion, i/)Tip, in our case is chosen to have 
vector quantum numbers so as to attach to an external 
photon. If we opt to use fermion bilinears for the meson 
interpolating fields, the possible Wick contractions fall 
into three classes: 

(a) the connected diagram of figure EJa); 

(b) disconnected diagrams where an intermediate state 
without charm quarks appears, e.g. figure Efb); 

(c) diagrams in which the photon couples to a closed 
fermion loop, e.g. figure Etc). 
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FIG. 5: Wick contractions for charmonium three-point func- 
tions, (a) connected, (b) OZI-suppressed, (c) disconnected 
loop 



Disconnected diagrams of the type in figure[S£b) should 
be small - in charmonium, OZI suppression is believed 
to be strong, which has at least a partial explanation 
in QCD perturbation theory since there the intermedi- 
ate two or three gluon state is suppressed as a power of 
a s {fx ~ m c ). We do not calculate diagrams of this type 
in the current computation. 

Diagrams with the photon coupling to a closed fermion 
loop, as in figure[5^c), are zero in the case of form-factors 
(i = /), on the basis of charge- conjugation invariance of 
the QCD action|2^|. In the case of transitions (i ^ /), 
while the "C-odd" connected two-point function vanishes 
in the ensemble average, the inclusion of the "C-odd" 
disconnected loop acts to compensate leaving a non-zero 
quantity. In the light quark sector such a diagram would 
be required to allow u> or <f> VMD-style processes, al- 
though such disconnected contributions to the QCD vec- 
tor correlator two-point functions are believed small as 
evinced by the small empirical a; — <f> mixing, and hence 
arc unlikely to contribute much to the transitions under 
study. In addition, in the particular case that the vector 
current is electromagnetic, we also have a suppression of 
this disconnected contribution owing to the sum of u, d 1 s- 
quark charges being zero and there being an approximate 
SU(3) flavour symmetry. In our quenched charmonium 
calculation we might consider allowing a charm-quark 
loop, but we expect the contribution to be small dynam- 
ically owing to the large charm quark mass and we do 
not explicitly include it in our computation. 

We are left with the connected diagrams, figure Ufa), 
with the possibility of light quark loops removed since we 
are working in the quenched approximation. We com- 
pute these diagrams on the lattice using the sequential 
source technique (2i|. This is simply a lattice technol- 
ogy that enables computation of the required propaga- 
tors with only two inversions of the Dirac matrix and 
without computing an all-to-all propagator. The con- 
nected three-point function written in terms of fermion 
propagators G° ! f(x,y) = i(jf '(x)ip? \y) is 

r /ri(P> 

= - E e-^e+^y tr(G(0, x)T f G(x, y)TG(y, 0) 
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In the sequential source technique, this is factorised as 

= _ J2 e+^y tr^Hl, (y, 0; t fl p f ) l5 T G(y, 0) r i75 ) , 
v 

with 

Hr f (y,0;t f ;p f ) = Y,^ rS G{y,x)^Y\ l5 G{x,Q) l5 . 



The cost is that we have specified the momentum and 
species of the sink particle and thus we must perform 
a new calculation for each new particle or momentum. 
For this reason we performed simulations only with pf = 
(000) and (100). We are able, however, in this formalism 
to insert any current, T, with any momentum, q. Our 
particular interest is with the vector insertions r = 7^. 
In addition, we can insert any meson source gamma ma- 
trix structure T; without further inversions. 



If the propagator G(y, 0) is computed via an inversion of 
the Dirac matrix, M^(x, z)GL(z, y) = 5ijS a p5 xy , then 
one can obtain Hy } {y, 0; tf\pj) via one further inversion 
using G(y, 0) in the inversion source, 

M(z, y)H Tf (y, 0; t f ;p f ) = 5 t „ tf e ip r%T\ l5 G(z, 0)75. 



The three-point function can be related to meson ma- 
trix elements of the vector current by inserting complete 
sets of states into iJSJ yielding 



,(3) 

m 



(p f ,q;tf,t) = e 2 Eftf f ) 7eM) (01^/^(0, 0)\f(p f ,r f ))(f(p f ,r f )\ijTij(Q, 0)\i(p it n)) ((0|^(0, 0)^^, n)))*, 



so that (in Euclidean space) we have a sum of decaying 
exponentials. If t is sufficiently far from and tf the ex- 
cited state contributions should have decayed away leav- 
ing us with only the ground state transition matrix ele- 
ment. The removal of the excited state contributions is 
assisted by using the interpolating field smearing choices 
we made to improve plateaux in the two-point functions. 

In our computation we place the source at t = and 
the sink at t = 24 and consider insertions of the vector 
current on all time-slices. Since we have an antiperiodic 
lattice of length 48, this makes the second half of the 
lattice symmetric with the first and in practice we av- 
erage the two halves to get our form-factor results, thus 
improving the statistics. 

The energies of meson states and the overlap of the 
interpolating fields with given states (Z) are extracted 
from the meson two-point functions as described in sec- 
tion IIIII This allows us to adopt a fitting method for 
extraction of matrix elements from the three-point func- 
tions. The procedure involves writing the three-point 
function on a given time-slice as the sum of products of a 
(known) propagation factor P = fj Z ^. e~ E f l f e~ < - Ei ~ Ef " >t , 
a (known) kinematic factor (e.g. = (pf +Pi)^ for the 
Vcy^Vc case) and an unknown form-factor, i.e. 

r(p/,Pi;i) = 'Y^P(pf,p i )t) ■ K n {p f ,pi) ■ f n (Q 2 )- 
n 

One can then form a vector of three-point functions that 
all have the same Q 2 = \pf — p, \ 2 — (Ef — Ei) 2 (i.e. ro- 
tationally equivalent momentum combinations and dif- 
ferent combinations of Lorcntz indices, labeled here by 



a, 6, c . . .) and a matrix with entries PK, to define a lin- 
ear system 



r(o;t) 




P{a;t)KM 


P(a;t)K 2 (a) 


T(b;t) 




P{b;t)Kr{b) 


P(b-t)K 2 {b) 


T( C ;t) 




P(c;t)ifi(c) 


P(c;t)K 2 (c) 



'h(Q 2 ) 



which we can invert with SVD to find the form-factors 

fn(Q 2 )[t}. ^ 

Single elimination jackknifc is implemented such that 
statistical errors are propagated through the inversion 
process. This proceeds by eliminating one configura- 
tion and averaging the remaining three-point function 
ensemble for use in the right-hand-side above and com- 
puting the propagation and kinematic factors using the 
equivalent "single-elimination bin" from the Z, E ensem- 
ble extracted from the two-point function fits. The SVD 
inversion is then performed for this "single-elimination 
bin" . This is repeated for elimination of each configura- 
tion yielding a configuration ensemble of solutions {/«}• 

The form-factors f n (Q 2 )[t] should have no t depen- 
dence once excited states have decayed away - i.e. they 
have plateaux. This linear inversion formalism can also 
include excited state contributions explicitly, provided 
their energies and Z's have been extracted from the two- 
point functions - we would simply enlarge the space of 
the linear system. In the present analysis we consider 
only ground state contributions - our poor estimation of 
the mass of excited states leads us to suspect that we 
are not able to isolate purely the first excited state. In 
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a few cases we will find that an f n (Q 2 )[t] does not reach 
a plateau within the available time, in such cases we fit 
in a region around t = 12 (away from the domain- wall 
oscillations near t = 0, 24), with a form 



fn(Q 2 



/n(Q 2 ) + f i e- mit + f/e- m ^ 24 -* ) , (9) 



with / ra (Q 2 ),fj,f/,mi,m/ as variables. A typical result 
of this fitting is shown in figure Efb). Our lack of control 
over the contribution of excited states is demonstrated 
by the extracted masses, uiij, which arc typically con- 
siderably larger than expected on the basis of either the 
true or our lattice first excited state masses. These fits 
are rather successful, suggesting that the single exponen- 
tial pollution hypothesis is reasonable even though we 
expect contributions from several excited states - a pos- 
sible explanation might arise if alternating excited state 
contributions oscillated in sign (cither from the Z-factors 
or from the excited state transition matrix elements) |3flj | . 
In the following sections, f(Q 2 ) points extracted in this 
manner are displayed in figures by ghosted shapes while 
points with a clear plateau appear with solid lines. 
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FIG. 6: Examples taken from J/ij) — > rj c J (a) plateau fits (b) 
fit by the function parameters are V{Q 2 ) = —1.55(1), 
fi = 1.45(3), fy = -0.42(14), ra, = 0.41(1), my = 0.27(7) 

A possible explanation for the large values of we 
observe comes from an alternate time-ordering allowed 



by our anti-periodic lattice. The periodicity implies that 
there is an image of the sink at t = —24, such that with 
the "current" insertion at timcslicc t > there is a pro- 
cess: creation of sink particle with Tf at t = —24, in- 
sertion of a "current" T,; at t = and annihilation of 
a vector particle by 7^ at t. This process has time de- 
pendence e~ 24E fe~ E *^ t , which would correspond in the 
ip — rj c case to = E^{q) ~ 0.5 in lattice units. A lim- 
ited computation of three-point functions with Dirichlet 
boundary conditions appears to confirm this hypothesis. 

In practice we can access a number of Q 2 values 
through projecting various momenta at the sink and the 
insertion. At the sink we are limited by the sequential 
source technique and only use pf = (000), (100), but at 
the insertion we project many momenta, corresponding 
to pi up to around (211). Higher momenta were calcu- 
lated but the signals were increasingly noisy. We used 
smearing at both source and sink with a = 3.6. 



A. Current renormalisation 

We utilise the simple local vector current, ^(x^^ip^x), 
which is not conserved with a discretised fermion action 
and as such gets multiplicatively renormalised by a fac- 
tor Zy{a). We extract Zy using a ratio of two-point and 
three-point functions evaluated at Q 2 = 0, which in the 
case of a temporal vector current is not polluted by ex- 
cited state transitions |3lT|. For the r\ c and the Xco we use 



±r (2) (n-f f 



24) 



where the factor of \ 



accounts for the equal contribution 
to the two-point function of the source at timeslice and 
the image of the source at timeslice 48 (recall that we 
have a temporally anti-periodic lattice). For the spin-1 
particles ( J/ip, Xci, h c ) we use 



4 M) (i) - 



IE, 



24) 



E(p)j: k r 



(3) 



(pf=Pi=p;t f = 24, t) 



Since we are working on an anisotropic lattice the renor- 
malisation of the spatial and temporal components of the 
vector current could, in principle, differ. The extracted 
Zy are shown in figure [7| where it is clear that there is 
no significant dependence upon the particle used in the 
extraction. The rj c extraction is cleanest and gives: 



Z ( y ] = 1.1803(2); Z. 



1.130(13). 



Note that these numbers differ by a factor comparable 
to the c values extracted in the dispersion relation anal- 
ysis, indicating that wherever a three-momentum p ap- 
pears we should in fact insert an explicit c = 1.032 to 
account for the imperfect dispersion. In what follows we 
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□ temp. p^(tOO) 
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Jib 



FIG. 7: extracted from two-point/three-point function 
ratios. Green dashed line shows Zy^ /c. The extracted spatial 
data are obviously compatible with this. 



will adopt this prescription for the momentum and use 
the temporal Zy. 

The renormalization factor, Za, for the four dimen- 
sional axial current operator = ipltils'tp is calcu- 
lated using the five dimensional conserved axial cur- 
rent for domain wall fermions A^ by the relation 
(A^A^O)) = Z A (A^t)A^0)). In the limit of L 5 -> 
oo, Zy = Za- However, at finite L§ there is a breaking 
of chiral symmetry that is characterized by a additive 
(and positive) quark mass shift called the residual mass 
TO ros . The difference of Zy from Za is proportional to 
this residual mass. The latter we find to be consistent 
with isotropic Wils on g auge calculations at the same spa- 
tial lattice spacing 



V. CHARMONIUM VECTOR FORM-FACTORS 

Charmonia, unlike charged or flavoured mesons, do not 
have radiative form-factors owing to them being eigen- 
states of charge-conjugation invariance. At the quark 
level this comes about through the photon coupling to 
both the quark and the anti-quark with equal strength. 
In our lattice computation we insert the vector current 
only on the quark line and are thus able to access its 
distribution within the meson. This can be compared to 
models. 



figure[S{b). Also shown is the expectation of VMD using 
the J/ip mass found in section UTTl this is seen to give a 
very poor description of the data in contrast to what is 
found in the light-quark sector for the pion form-factor 
p9j . This is to be expected on the basis of a disper- 
sion equation approach - in the light quark sector, the p 
pole is the nearest left-hand vector singularity in the com- 
plex Q 2 plane and the next nearest (neglecting multi-pion 
cuts and isospin suppressed poles) is the excited p'(1450) 
which is relatively distant. In charmonium the nearest 
pole is the J/V'(3097) and the next nearest the ^'(3686) 
which is not sufficiently distant as to be negligible. Hence 
in charmonium it is likely that one needs to sum many 
vector meson poles (of unknown residue) to agree with 
Q 2 > data. 




'b 
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FIG. 8: 7y c form-factor (a) typical lattice signal with plateau 
fits (b) form-factor 



A. rj c "form-factor" 

This is defined by the matrix element decomposition 

yic{p f )\f{Q)\vM)) = f(Q 2 )(p*+P f r- (10) 

In figure 13 a) we display some typical f(Q 2 )[t] as ex- 
tracted from our simulation in the manner described in 
the previous section. We see clear plateau behaviour 
away from the source and sink points which is fitted and 
a value of f(Q 2 ) obtained. These points are displayed in 



We take the more pragmatic approach of fitting 
the (Q 2 > 0) data to a simple analytic form, 

cxp — -Op? (l + ctQ 2 ) ■ This fit is also shown in figure 

ISJb) where one can see it faithfully reproduces the data 
over the entire Q 2 range considered. The fit parameters 
so obtained are 

[3 = 480(3)MeV; a = -0.046(l)GcV~ 2 . 

One can define a squared "charge" radius by (r 2 ) = 
-6^r/(Q 2 )| Q2=0 = which yields ^/J^ 2 ) = 
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0.255(2)fm, giving some a posteriori justification for our 
lattice size being only ~ 1.2fm. 

As mentioned in the introduction, most charmonimn 
phenomenology is done within the framework of quark- 
potential models. These models are usually of the non- 
relativistic Schrodinger equation type, utilising some 
phenomcnologically justified static potential and incor- 
porating relativistic corrections to this, such as spin-spin 
and spin-orbit terms. By and large these models are suc- 
cessful so it makes sense for us to compare our lattice 
results with them in cases where there is not experimen- 
tal data. 

In the simplest forms of these models we have only 
Galilean invariance and not the full Lorentz invariance 
of a relativistic theory. This is manifested in the factori- 
sation of the spatial wavefunction of a charmonium state 
into centre-of-mass and internal wavcfunctions: 

\ Vc (p)) oc / d 3 Re^ [ dVV(r), 



[ kT" 




f ° 2 1 


16/3 2 


= exp 


16/3 2 _ 



with ip(r) independent of the r\ c momentum p, obtained 
by solving a Schrodinger equation in the rest frame of the 
rj c . This poses problems of frame-dependence when one 
calculates matrix elements. Consider calculating f(Q 2 ) 
using equation (|1L)|> within this model. We want to work 
in a frame near to the one in which the r\ c wavefunction 
was calculated, to minimise the neglected boost effects, 
two choices strike us: 

• the rest frame of the initial r\ c (pi = 0, pf = q). 
Here only the end-state rj c could suffer boost dis- 
tortion. We have 

/(Q 2 = 2m r 2 J-l+(l+^) 1 /2]) oc J d^r(r)j^)m 

With approximate harmonic oscillator wavcfunc- 
tions (tp(r) oc e _/3 r ) we obtain 



of the matrix element will be essentially identical 
to the previous case with the only difference being 
that now Q 2 = \q\ 2 . Hence 



/bf(<2 2 ) -> exp 



and we explicitly see the frame-dependence of the 
simplest quark models (/bf(Q 2 ) ^ /rf(Q 2 ))- 

If Q 2 is small, all frames will agree, and it is reall y on ly 
here that we have a unique quark model form. In [33 , a 
thorough analysis of these concerns is applied to the case 
of the proton form- factors. 

We note in passing that light-front quark models have 
the same sort of problems - although boosts along a given 
axis become simple (kinematical), which may admit con- 
sistency up to higher Q 2 , we lose rotational invariance, 
and this can lead to extra form-factors not present in a 
picture with 0(3) symmetry |34j. 

In the quark model of ISGW [35|, form-factors have 
a dependence near Q 2 = of the form exp [ — jgjp^r ] , 
with k ~ 0.7. The k factor was added by hand to better 
describe the pion form-factor and certain heavy quark 
transitions and was ascribed an origin in relativistic cor- 
rections 5 . An alternative origin might be the effect of 
gluonic degrees-of-freedom not considered in the simple 
potential model jU H3- In a Schrodinger equation 
with a Coulomb plus linear potential and spin-dependent 
corrections was solved using a variational harmonic os- 
cillator basis - they found for the r/ c that (3 = 710McV. 
Hence k(3 = 500MeV which is in rather good agreement 
with our (3 = 480(3)McV, suggesting that the potential 
■model, with the relativistic correction, is capturing at 
least some of the Q 2 « physics we have on the lattice. 



/rf(<2 2 ) ->■ exp 



f kT" 




16/3 2 _ 


= exp 



16/3 2 



Here we share 



• the Breit frame (pi = — |, pf = |) 
the possible boost distortion equally between the 
initial and final states. The quark-model evaluation 



B. J/ip "form- factors" 

A vector particle has three form-factors, a suitable set 
being those of definite multipolarity: charge - Gc(Q 2 ), 
magnetic dipole - Gm(Q 2 ) and quadrupole - Gq(Q 2 ). 
These are defined in terms of the standard Lorentz co- 
variant decomposition 



(V(p f ,r f )\f(0)\V(p u n)) = -(p f +p i r[G 1 (Q 2 )e*(p fl r f ) ■ e(p u n) + 201? $ fl n ). Pl e(p h n).p f 

+ G 2 (Q 2 ) e ll (p l ,r i )e*{p s ,r f ).p l + e^{p f ,r f )e{p l ,r i ).p s 
I 



by 



G c = (1 + ft;) Gi - §77G 2 + §77(1 + v)G 3 
Gm = — G2 

G Q = Gx - G 2 + (1 + 7?)G 3 , 



O 1751 

with rj = 2^-2 . In figure |£] we display our lattice points 
and a fit to them. The charge multipole is fitted by the 
same function we used for the r\ c with resulting paramc- 
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FIG. 9: J/tp "charge" form-factor 



tcrs: 



pc = 470(7)MeV; a = -0.022(7)GeV" 



(11) 



This similarity of these parameters to those extracted for 
the r/ c indicate that these two particles have spatial wave- 
functions that are rather alike. The magnetic dipole and 
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FIG. 10: J/ip form- factors (a) magnetic (b) quadrupolc 
quadrupole form-factors, displayed in figure ITU1 have fits 



of the form G(0) exp - 
ful and yield 



16/3- 



which are seen to be success- 



M 



520(8)MeV; 
580(44)MeV; 



G M (0) -- 
Gq(0) 



2.10(3) 
= -0.23(2). 



The quadrupole form-factor at Q 2 = gives us the 
quadrupolc moment of the J/tp and, via a quark- 
model interpretation, access to the quark-antiquark D- 
wave admixture in the dominantly S- wave J/tp. The 



quadrupole moment operator is [4fJ Q = r (3 cos 9 — 

1) = J^r 2 Y§{6,(p). Then it is clear that (S\Q\S) = 0, 

while (S|Q|Z?) 7^ 0. In the approximate harmonic oscil- 
lator basis, we have 

(J/V>|Q|JA/>) = 2V3a s a D P~ 2 , 

where is the L-wave amplitude in the J/tp wavefunc- 
tion. Since Gq(0) = m 2 ^(Q}, using our lattice J/tp 

mass we have (Q) = -0.027(3)GeV" 2 . With the ex- 
tracted quadrupole (3q value, and assuming as ~ 1 we 
find cld = —2.6(5) x 10~ 3 . This very small value can be 
used to constrain the size of any tensor term one might 
wish to add to a charmonium model Hamiltonian. We 
remind the reader that any contribution from coupled- 
channel DD loops will not be present in our quenched 
computation. 

The magnetic dipole moment, /ij/^, can be extracted 
using Gm(0) = 2Mj/4> l^j/ip - hi units of the J/tp mass 
it is 1.05(2). Within the simple non-relativistic quark 
model this quantity can be expressed as 

(1 + 5) [(1 + K c )\a s \ 2 + |(1 - 2 Kc )\a D \ 2 ] , 



where S 



M, 



-2m c 



and k c is the anomalous magnetic 



moment of the charm quark. With the D-wsve admix- 
ture extracted from the quadrupole moment, the second 
term in square brackets is negligible. If one had a way of 
unambiguously setting the charm quark mass one could 
use this expression to extract the charm quark anomalous 
magnetic moment - however this is not possible even in 
principle - the quark mass is a renormalisation scheme de- 
pendent quantity. In quark potential models it is usually 
tuned along with other parameters to give a good de- 
scription of the spectrum. A conservative interpretation 
of our extracted value would be that it could be explained 
in a model with a small or zero anomalous quark moment 
and a charm quark "mass" slightly less than half the J/tp 
mass. This would be interesting application of the EFT 
method of pNRQCD [41|- within this picture one has, at 
a given order in the power counting, a relationship be- 
tween the charm quark mass and the J/tp mass in terms 
only of the (determined) strong coupling, as such one can 
determine, at the same order in the power counting, the 
anomalous moment k c . 
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C. \co "form- factor" 

This has the same decomposition as does the rj c , eqn 
1)10(1 . Our lattice three-point functions here are much 
noisier than for the r\ c and J/ip, but there are sufficiently 
clear plateaux to extract the values plotted in figure ITT1 

_ Q 2 

The data is fitted with the form /(0)e 16 0' 2 giving /(0) = 
1.0015(15), consistent with 1 as it should be, and (5 — 
393(12)MeV, which is roughly consistent with the quark 
model value |39| k(3 = 340MeV. That this is smaller than 
the r) c , J/ip (3 is simply a reflection of the fact that the Xco 
is spatially larger due to the centripetal P-wave barrier 
between the quark and the anti-quark. 




FIG. 11: XcO form- factor 



VI. RADIATIVE TRANSITIONS 

A. J/tp — > rjcjMi 

The Minkowski space-time matrix element for this 
transition can be expressed in terms of a product of one 
Lorcntz invariant form-factor and one Lorentz covariant 
kinematic factor: 

(Vc(p')\nO)m(P,r)) = 2Vi f ] e^p' aPp e,(p,r). 

(12) 

This decomposition is parity invariant and, with V(Q 2 ) 
real, time-reversal invariant. The Lorentz invariant ma- 
trix element for J/t/j — > ?7 C 7*(Q 2 ) is the contraction of 
(11211 with a final-state photon polarization vector: 

M r ,r y = e^r^^WmJ/mr))- 
The decay width with an on-shcll photon is 



q 32ir 2 ml 3 



(13) 



which contains a sum over the final state photon polar- 
isation and an average over the initial ip polarisation. 
Explicitly then we have 



T(J/ip -> ?7 C 7) 



1 



kTl 



47r (m 



In our lattice computation we couple only to the quark 
and not to the anti-quark and do not include the quark 
electric charge factor. As such we compute V(Q 2 ) which 
is related to V(Q 2 ) by 6 



V(Q 2 ) = 2x -exV(Q 2 ), 



and hence 



T(J/ip -> i] c j) = a 



{m Vc + m^) 2 27 



6 V(0)| 2 . 



(14) 



There is only one direct experimental measurement of 
this width [H, r CB (J/V' -> Vcl) = 1.14(33)keV. There 
arc, in addition, several measurements of the product 
branching ratio for the process J/ip — > rj c j — > _^rf>70 and 
one independent measurement of rj c — > cfxj) |43|. Taken 
together these data give Y^{J/\jj — > r] c j) = 2.9(1.5)kcV, 
which is consistent within the large errors with the Crys- 
tal Ball result. 

There is an essential ambiguity in how we compare our 
lattice results with the experimental data which arises 
from the fact that our lattice charmonia masses do not 
coincide exactly with the physical masses. The prob- 
lem lies in whether we should use experimental or lattice 
masses in equation (|14fl. |g| is closely related to the hy- 
pcrfinc splitting which is rather sensitive to details of the 
lattice calculation (as discussed in section InTf . hence we 
observe considerable difference in using lattice or exper- 
imental masses: 



mo)i 



expt. 



Crystal Ball 



phys. masses 1.27(19) 2.02(52) 
latt. masses 1.56(22) 2.48(64) 

We extract V{Q 2 ) from three-point functions comput- 
ing using the following sequential sources at the sink: 
tt[p) = (0,0,0); (l,0,0)];p ( x, tf ,,)|?/ = (0, 0, 0); (1, 0, 0)]. 
The values from different sink particle and momentum 
choices are all seen to be consistent suggesting that the 
sequential source technique is working correctly. The lat- 
tice data displayed in figure IT21 is fitted with 



V{Q 2 ) = V(0)e 



Q z 



where we find 



V(0) = 1.85(4); /3 = 540(10)MeV. 



6 V(Q 2 ) is what would appear if we were computing the 
transition in the isospin limit 



' 7T 1 7 
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This agrees reasonably with the Crystal Ball result in 
the case of using lattice masses throughout. Without 
performing a computation at the correct quark mass and 
relaxing approximations sufficiently to duplicate the ex- 
perimental hyperfine splitting, we cannot give a more 
definitive result. 



At small Q 2 , \ff - {m \^ ? + Q 2 \(l + $A , which 

means the simple quark model with harmonic oscillator 
wavcfunctions has the same Q 2 -dependence as our fit. 
Other reasonable frames, in the small Q 2 limit, have the 
same dependence. 

Performing the full calculation (with k c = ) at Q 2 = 
in a quark model gives an expression like |44l |45| 
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FIG. 12: J/tp — » ry c 7 transition form-factor 

The gaussian fit performed to the data was clearly 
rather successful; if we wish we can motivate such a form 
by returning to the sort of simple quark models consid- 
ered in the previous section. Within these models one 
performs a non-relativistic reduction of the current op- 
erator and computes the matrix element of this between 
the J r /tp and r\ c wavefunctions. This is a magnetic dipole 
(Mi) transition that occurs through quark spin- flip - the 
lowest order current operator (in an expansion in v/c) has 
the form c?x {p' q —p q ). Computing in any frame (recalling 
that the wavcfunctions are assumed to be unchanged un- 
der boosts and that we have the minimal non-relativistic 
reduction of the current) we find 



3m 2 



r 2 drR; c (r)j (^)R^r 



which can be used to evaluate the width provided the 
wavefunctions and the mass of the charm quark are 
known. Determining this latter quantity is a tricky is- 
sue - it is usually set at the Schrodinger equation stage 
along with the phenomcnological potential parameters 
such that the mass spectrum is reproduced, however a 
unique solution is not easily found and as such this mass 
parameter can vary considerably between mod els, e.g. 
m c = 1.84GeV in Q while m c = 1.628GeV in and 
even m c = 1.479GeV in |4^. Given this uncertainty, es- 
timates for the Mi transitions are rather imprecise, de- 
pending as they do, on the inverse square of the mass 7 . 
In principle we do not have this problem in our lattice 
simulations - while we do vary the quark mass to get 
agreement with the spectrum (we have to set the charm 
quark mass somehow) , we have no other free parameters 
to vary. Unfortunately, as discussed in section UTTl we did 
not tune the quark mass perfectly and this is a source 
of some systematic error on our result for this transition 
and may be the reason our result is somewhat too large 
(recall that our quark mass is slightly too small). 



V(Q 2 )<x / r 2 drRl(r)j (^)R4r) 



exp 



16/3 2 



B. 



XcO 



for harmonic oscillator wavefunctions. If the rj c and J/ip 
are allowed different f3 then (3 2 = \((3 2 a + Once 
again we encounter the problem that V(Q 2 ) should be 
frame independent, but \q\ is not. In the rest frame of 
a decaying J/tp, a "natural" frame to consider for this 
process, 



(m 2 



l lf + 2Q 2 {m%+ml) + 



4m 2 



Experimentally it is only possible to access transverse 
on-shell (Q 2 = 0) photons in this transition and the ma- 
trix element is purely through the electric dipole {E\ ) . In 
more generality, if we allow Q 2 ^ and necessarily also 
longitudinal photons, there is a second multipole, labeled 
C\. The decomposition of the transition matrix element 
in terms of these multipoles is derived in the appendix 
and we reproduce it here: 
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(S(psW(0)\V(Pv,r)) = IT^Q 2 )! £ x (Q a )[n(QV(fr.O - e(p v ,r).p s (p»p v .p s - m 2 v p») 

, Ci(Q 2 ) 



-m v e{p v ,r).p s 



Pv-psipv +psY ~ "isPy - m vPs 



The Lorentz invariant matrix elements for the transition 
J I 'ipj* (Q 2 ) are also given in the appendix: 



XcO 



T)=E 1 (Q 2 ) 
0;r> = 0) = -Ci(Q 2 ) 



M{r 1 = 
M(r 7 

Hence the analogue of (|13|l gives for the width at Q 2 = 0, 



a- 



\q\ 



^(o)| 2 , 



XcO 



where the lattice form-factor is again related to the phys- 
ical one by E X (Q 2 ) = 2 x |e x £i(Q 2 ). 

The most recent measurement of this decay's branch- 
ing fraction comes from the CLEO collaboration^, who 
find, using the PDG total width to normalise: r(xc-o ~ ¥ 
J/ipj) = 204(31)keV. In addition to this we have the 
PDG|3] average/fit to data obtained up to 2005 which 
gives r(x c o -> J/th) = H5(14)keV. The next PDG re- 
port will likely contain the CLEO value in a new average 
which will thus lie between these two values. 

In figure IT31 we display the Ei(Q 2 ) extracted from our 
lattice simulations. Temporal vector current insertions 
produce compatible results but with much larger error 
bars and are not shown. 

Our simulation data lies at Q 2 ^ 0, but since we are 
primarily interested in the photopoint we require some 
fit function to allow us to extrapolate back. In the light 
of the success of forms motivated by the non-relativistic 
quark model in previous sections we consider using a 
function which resembles one that would be derived in 
such a model. We opt to use a form 



E 1 (Q 2 )=E 1 (0)[1 



cxp 



Q 



2 1 



16/3 2 



(15) 



which has the gaussian behaviour used previously modi- 
fied by a polynomial in Q 2 . In the simple quark model, 
the Q 2 / p 2 term could arise from relativistic corrections or 
departures from gaussian wavefunction behaviour. Note 
that this form is analytic for Q 2 > as we would expect 

- singularities (as in the VMD case) will occur at Q 2 < 0. 

We do not include in the fit the points at Q 2 < 

- these data, corresponding to the case pf — pi where 
Q 2 = —(Ef — Ei) 2 , were extracted from correlators with 
no plateau behaviour using the fitting method described 
in section IIVI It is therefore a rather non-trivial cross- 
check that our fit function, constrained by points at Q 2 > 
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FIG. 13: XcO — * J/ipl Ei transition form-factor, (a) full range 
of lattice data (b) zoom to the Q 2 ~ region 



lGeV 2 , extrapolated to the Q 2 < region, overlays these 
points. 

The fit returns the following parameters: 

a t £i(0) = -0.137(12) 
P = 542(35)MeV; p = 1.08(13)GeV 

The longitudinal photon transition form-factor, 
Ci(Q 2 ) can also be extracted from lattice three-point 
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functions. We display our data in figure [21 As dis- 
cussed in the appendix, Cl ^_ - 1 is real for real Q 2 = — q 2 if 

time-reversal invariance holds - thus for Q 2 > 0, C\(Q 2 ) 

_Qf_ 

is imaginary. The fit is to the function icy Q 2 e 16 ^ 2 
which has the required property that Ci(Q 2 — > 0) — > 
and which is what one would expect in a simple quark 
model. 



The fit suffers from the large error bars on the data, 
but does at least yield a j3 value compatible with the 
value extracted from the E\ fit: 

a t c= 1.83(16); (3 = 501(33)MeV 
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FIG. 14: XcO — * J/ip^f Ci transition form-factor 



With real photons this transition receives contributions 
from two multipolcs, the dominant electric dipole (Ei) 
and a much suppressed magnetic quadrupole (M 2 ). Ex- 
perimentally the M 2 contribution is measured through 
the angular distribution of photons - the PDG[]| aver- 
age the two extractions performed^i^^i each of which 
found a number consistent with zero, to give 



M a (0) 



V^i(0) 2 + Af 2 (0) 2 



-0.002jft8» 



Appendix ^ contains the tools required to derive the 
relation connecting the transition matrix element with 
the multipole amplitudes: 



(A(f A ,r A )\f(0)\V(p v ,r v )) 



4V20(Q 2 



'(PA -Pv)a* 



E 



(Q 2 )(PA + pv) P {^m A [e* (pA,r A ).pv] ^ (Pv ,r v ) + 2m v [e(p v , r v ) .p A ]el (p A , r A ) J 
+ ^/2(<3 2 )(PA +py) p ^2?Ti J 4[e*(pA,?'A)-Py]e 1 y(pv,7'v) - 2my [e(pV, ry).p yl ]e* (px, r A fj 



C1(Q 2 



- ^VL(Q 2 )el(p A ,r A )e p (p v ,r v ) 



+ (pa+Pv) p {m A -my + q 2 )[e*(p A ,r A ).p v ] e v (p v ,r v ) + {m 2 A - m v - q 2 )[e(p Vl r v ).p A ] e*(pl,r A ) ) 



(16) 



In extracting these amplitudes from our lattice three- 
point functions we are struck with the problem that our 
Xd signal (from the operator 7*7 5 ) becomes noisy after 
relatively few timcslices, in fact at roughly the same time 
that the two-point function begins to plateau (see figure 
[Hb,c)). Because of this it was only possible to extract 
a convincing signal for the Xcl with Pf — (000), which 
considerably limits the number of three-point functions 
available to us for the transition Xcl — * J/ip"/- In ad- 



dition the plateaux in the few available correlators are 
borderline, so that for the lowest Q 2 point a fit of the 
form © was applied. 

In figure [TBI we show the extracted E\(Q 2 ) and a fit of 
the same form as used for Xco ~ * J/ipjEi ■ Within the 
large extrapolated error we are in agreement with the 
experimental data. Unfortunately, unlike in the Xco 
J/tpJ case, the E\ kincmatical factor for pf = pi where 
Q 2 f» 0, is zero, so we cannot directly obtain E\(Q 2 w 
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0) and the subsequent cross-check on the extrapolation 
which gave parameters, 

Ot^i(O) = 0.312(39) 
(3 = 555(113)MeV; p = 1.65(59)GeV 
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FIG. 15: Xci - * J/ipJ Ei transition form-factor 

In figurclTrjlwe show the extracted M^iQ 2 )- The fit has 
the same functional form as in the E\ case and returns 
parameters: 

a t M 2 (0) = -0.062(37) 
f3 = 617(142)MeV; p = 0.93(47)GeV 

The (3 value is compatible within the large error with the 
value extracted from the E\ data. 




2 3 

d (GeV 2 ) 



We observe in figure HBl that our lattice data extrapo- 
lates down to a value consistent with zero and hence with 
experiment, within large errors, at Q 2 = 0. 

There are a number of approaches we might consider 
to reduce the error bar on the predicted value. The 
first is simply brute force; evaluate on a larger num- 
ber of gauge field configurations, thus reducing the sta- 
tistical fluctuations. Another option, which would also 
require increased computation time, is to work with a 
larger spatial volume (at the same lattice spacing); this 
allows smaller discrete three-momenta and hence access 
to points closer to Q 2 = 0, the reduced extrapolation dis- 
tance then reducing the error on the extrapolated point. 
A third possibility, ideally combined with the previous 
two, is to enhance the Xd plateaux by finding an in- 
terpolating field with maximal overlap on to the ground 
state Xci- This might involve diagonalising a matrix of 
two-point functions in a basis of different smearings. 

Finally we have the longitudinal photon multipolc C\ 
which we display in figure 1171 along with a fit of same 
type as used for the Xco — * J /''Pi transition. From the fit 
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FIG. 17: Xci — * J/ipf Ci transition form-factor 
we extract 

o f c= -2.91(26); = 502(38)MeV 



D. h c 



FIG. 16: Xci - * J/V'7 M2 transition form-factor 



The h c , with J PC = l" 1 , was only recently ob- 
served with high significance by the CLEO collaboration 
[49I IHflj . The reason for the delay of discovery with re- 
spect to the other ground state charmonia lies in the dif- 
ficulty of production - the process eventually utilised was 
the isospin violating ip(2S) —> ir°h c with a subsequent 
h c — > rjcl- Since only the product branching fraction 
B(tp(2S) — ► TT°h c )B(h c — > r) c -f) is measured, our calcula- 
tion of T(h c — > Vcj) constitutes a prediction. 
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The form-factor decomposition is identical to the 
XcO ~ *■ J/^l case. Here again we suffer from poor two- 
point functions - we were only able to extract a convinc- 
ing h c signal in the Pf = (000) case and as such we can 
use only three-point functions with the h c at rest. De- 
spite this limitation we are able to extract some clean 
plateaux in a few cases and in the remaining cases we fit 
using 10 - we plot the extracted E% in figure^] The fit 
shown (which is to all points including the one at Q 2 < 0) 
is of the form 



£i(0)exp- 



16/3 2 



since the data does not seem to require a term (p 1 
in such a fit is very small). We find 



-0.1 



-0.2 
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FIG. 19: h c — > rjc'y C\ transition form-factor 
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FIG. 18: ft c — > ??c7 Si transition form-factor 
The relation to the width is 

r(/j c ^ c7 ) = a#^i(0)| 2 , 

so that we predict 

r (/lc ^, c7) = 663 ( 132 W 

K 1 " 601(55) 

where the upper value uses lattice masses and the lower 
value physical masses. 

The Ci multipolc is displayed in figure UU with the 
same fit form used for the XcO ~ > J 1^1 transition. 



a t c 



-2.90(18); /3 = 545(49)MeV 



VII. DISCUSSION 



We extracted, for the first time, a (limited) char- 
monium spectrum using domain wall fermions on an 
anisotropic lattice. While an imperfect quark mass tun- 
ing left our spectrum systematically 5% too light, the 
gross features of the S & P levels were correct and no- 
tably we found a rather large hyperfine splitting in con- 
trast to other quenched improved actions. 

By attaching a vector current to only the quark line 
and not the antiquark, we avoided the constraint of 
charge conjugation invariance and sampled the vector 
form-factors of the lightest three charmonium states. The 
charge form-factors of the rj c and J/tjj were similar as 
expected, and, as one might anticipate on the basis of 
dispersion relations, did not appear to be dominated by 
a single VMD mechanism. 

For the J/ip the magnetic dipolc form-factor was ex- 
trapolated back to Q 2 — to yield a magnetic dipole mo- 
ment that was consistent with there being zero anoma- 
lous charm quark magnetic moment in a quark model 
picture. Similarly the very small quadrupole moment 
extracted indicated minimal D-wave admixture into the 
J/tp- The XcO form- factor displayed a larger charge ra- 
dius than the r/ c , indicating the effect of the centrifugal 
barrier in the P-wave state. 

In Table [H] we summarise our results for radiative 
transitions, where it is clear that the (3 values for the 
P —> S transitions are all compatible with having the 
same value, which is in line with quark model expecta- 
tions that the Xcj,h c spatial wavefunctions should dif- 
fer only through small spin-orbit distortions (and per- 
haps differing closed-channel effects if these are allowed 
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for[s3 |53 8 ). Within (p)NRQCD one might also have 
differing color octet contributions |4lj . 

El XcO — ► J/ipf Xd — > J/ipl h c — » r; c 7 

P/MeV 542(35) 555(113) 689(133) 



p/MeV 1080(130) 1650(590) 



p lat.mass /UpV 



■ phys.mass 

-» PDG 
CLEO 



rPRG/keV 



288(60) 
232(41) 
115(14) 



600(178) 
487(122) 
303(44) 



OO 
663(132) 
601(55) 



Ml 


J/l/>-> Vcl 


M2 


Xd -* J/</>7 


/3/MeV 


540(10) 


/3/MeV 


617(142) 


p lat.mass ;^.y 
phys.mass/ v 


1.61(7) 
2.57(11) 


M2 
El 


-0.199(121) 


r™ G /keV 


1.14(33) 
2.9(1.5) 


expt. 


-0.002(+ 1 8 7 ) 



CI XcO — > J/V>7 Xd — > J'/V'T — * ??c7 
/3/MeV 501(33) 502(38) 545(49) 
|£|/GeV 11(1) 17.6(1.6) 17.5(1.1) 

TABLE II: Radiative transitions 

We can also compare the pattern of p values for the 
Ei transitions with the expectations of a simple quark 
model. Performing a non-relativistic reduction of the 
vector current, the p term arises from the spin-dependent 
correction (oc a x q) to the dominant convection current 
oc p; using effective harmonic oscillator wavefunctions 
one finds (531 



^ M (Q 2 ) = a|l + ^,exp. w ., 



(17) 



where r is related to spin-orbit Clebsch-Gordan coeffi- 
cients, 



errors on the lattice results, these relations appear to be 
satisfied. 

Within the quark model, the transition is sup- 
pressed relative to Ei by one power of v/c. It is also 
rather sensitive to any charm quark anomalous magnetic 
moment. Some details are worked out in |54| . where they 

find a value (setting k c = 0) M 2 (0)/J5i(0) 0.06. Our 

data is unfortunately not sufficiently accurate to discrim- 
inate on this level - we outlined earlier in the text some 
possible improvements to the calculation to remedy this. 

In this first attempt at charmonium radiative transi- 
tions using lattice QCD we have demonstrated that it 
is possible to get reasonable agreement with experiment 
and have gone some way to justifying certain results of 
the more widely applied quark model. Future lattice 
work in this direction will have to address the problem of 
reliable excited state extraction in order to consider such 
well-measured transitions as ip' — > XcJl- 

There is, naturally, a desire to see calculations done 
without the quenched approximation, but, as discussed 
in section [HJ we do not expect unquenching to affect 
radiative transitions particularly strongly, except in the 
sense that it will improve the lattice state masses and 
help remove the phase-space ambiguity we encountered 
in section IVII However, an unquenched computation is 
warranted to test models which propose a considerable 
effect from coupled channels |5lj. 

Our ultimate aim is to study photocouplings of light- 
quark hybrid mesons, with this in mind the next step 
will be to consider radiative transitions involving char- 
monium hybrids - the non-local interpolating fields re- 
quired for this study will also allow us to access higher 
spin conventional charmonia such as the \c2- 



2 XcO 
r = 1 Xd- 

h e 

Working in the Xc rest frame at small Q 2 we would have 



1^ « k1^ + Q 2 (l + A) where \g\ - ^ 

momentum transfer at Q 2 = and A = i ^ 
Thus we can express the quark model form as 



is the three 
- 1 



E? M {Q 2 ) = E? M (0) 1 



1 + A 



4/33, 1 



exp - 



with S = r|]f and ^ M (0) = a{\ - 



-*)exp-Jgk. 



Q 2 (l + A) 
16/3 2 : 

(18) 
Hence, 



so 



to a first approximation we'd expect that p 

that p(xci) ~ V%p(Xco) and p(h c ) — ► oo. In the same 
approximation we have p(Xco) ~ 2/3^. Within the large 



But note that such effects are not present within our quenched 
calculation 
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APPENDIX A: MULTIPOLE DECOMPOSITION 

It is convenient to express radiative transition ampli- 
tudes in terms of multipoles. In this appendix we derive 
Lorentz covariant decompositions of vector current ma- 
trix elements into multipoles. These decompositions do 
not appear to have been explicitly presented previously 
in the literature. 

Our method involves writing down the most general 
Lorentz covariant, current conserving and parity invari- 
ant decomposition of the matrix clement of the current 
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in terms of a number of arbitrary form-factors. We then 
compute the helicity amplitudes for the decay i — > f"f 
by contracting the current matrix clement with a photon 
polarisation vector. For convenience we work in a partic- 
ular frame, but the result is covariant. The relationship 
between helicity amplitudes and multipolcs is prescribed 



in |55| , whence we eliminate the arbitrary form-factors in 
favour of the multipolc form-factors. 

We will demonstrate the method with the scalar- vector 
transition (0 + «-> 1~) relevant to Xco - * J/^l- The most 
general Lorentz covariant decomposition is: 



(S(psW(0)\V(pv,r)) 

= A(q 2 )e^e p (p v ,r)p Vr7 p ST + P(q 2 )e»{p v ,r) + [e(p v ,r).p 3 ] [B + (q 2 )(p s + p v Y + B_(q 2 )(p s - p v Y 

(Al) 

where the polarisation vectors carry a label r which is the z-component of the spin (which is not equal to the helicity 
in general). 

Parity invariance requires that 

(S(psW(0)\V(p v ,r)) = (Sips^V-'Vj^O^VlVipv^)) 

= -PZ(S(-ps)\r(0)\V(~p v ,r)), 

where we've used T\S{p s )) = \S{-p s )), V\V{p v ,r)) = -\V{-p v ,r)) and V^j^V = P£?" 9 . Hence 

A(q 2 )e^ aT e p (p Vl r)p Va p S r + P(q 2 V(pv,r) + [e(p v ,r).p s }[B + (q 2 )(p s + PvY + B-(q 2 )(p S - PvY 
= -P^A(q 2 )e^ T e p ^p v ,T){Pp v ) a {Pp s ) T + P(q 2 )^{~p v ,r) 

+ [e(-p v> r).(Pp s )] [B + (q 2 )(P Ps + Ppyf + B_(q 2 )(P Ps - Pp v Y 

I 



Properties of the rotation group give that e A '(— p, r) = 
— P£e v (p, r) using which one verifies the parity invariance 
of this decomposition provided A{q 2 ) = 0. 

The conservation of the vector current is an additional 
constraint, 

Q = d^S(p s )\f(x)\V(pv 7 r)) 

= d^S(p s )\e^ x r(0)e-^\V(pv,r)) 

= d^ Ps - p ^- x {S(psW(0)\V(f v ,r)) 
= q p (S(p s )\f(0)\V(pv,r)), 

whence we eliminate one of the three form-factors, 

P(q 2 ) = (m 2 v - m 2 )B + (q 2 ) - q 2 B_(q 2 ). 

Multipolc amplitudes are most easily defined in terms 
of helicity amplitudes which can be obtained from the de- 
compositions of the previous section in a straightforward 
way. First we will find the defining relation for multipolc 
amplitudes in terms of helicity amplitudes. This is done 
in analogy with Durand [55| , but rather than working in 
the Breit frame we choose to work in the rest frame of 
the decaying particle. The results are Lorentz covariant 
so this choice of frame is irrelevant. 

Wc can define a vertex function in the rest frame of a 



particle of spin- J 

r"(J'A'; JA) = (p'z;J'X'\f\JX) 

= (J'\'\e^' K3 f\J\), (A2) 

where (p'z\ J'X'\ is a final state of spin- J', helicity-X' 
in motion along the positive z-axis with momentum p' , 
| JA) is a similar (initial) state at rest and e l ^p ,I<3 is the 
unitary operator effecting the boost from rest to mo- 
mentum p' z. The matrix element <|Alll discussed in the 
previous section is a generalisation of this vertex func- 
tions to an arbitrary Lorentz frame, i.e. r M (0,0; 1, A) = 
(S , (p'z)|j' i (0)|y(0,r = A)) . Note that now we are opt- 
ing to use helicity and not the z-component of angular 
momentum. However since we have chosen the end-state 
particle to move along the positive z-direction, we have 
A = r, A' = r'. 

In this frame the amplitude for the decay (JA) — > 
(J' A') + ( 7 A 7 ) is 

M = e* v (—p'z, A 7 )r^( J'A'; JA), 

where the polarisation vectors for a photon moving along 
the (— z)-direction are 

e< l (-p'z,X y =±) = ±j=(0,l,Ti,0) 

e^-p'zA, = 0) = -4- 
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We find that 

M(X 7 ) = <SA,v-A 7 CA,r^(J'A'; jx), 

# and T± 



where c± = l,Co = and T± = T^r 1 ± iT 2 ). 

Returning to the definition of the vertex function l|A2|l . 
we assert that e**>> s j I/ can be expressed as a sum of 
operators which transform as tensors under the rota- 
tion group - this is the essence of a multipole expansion: 
e % *p' K3 j v = T^ v y The vertex operator is then a sum 
of matrix elements of these tensors, each of which sat- 
isfies the Wigner-Eckart theorem, so that, for all spins 
integral, we can write 



v(j'\';j\) = Y / (-iy'+x ^ * 



{J'\K } \\J), 



where v = ±, 0. j° and j -± are in different representations 
of the rotation group, so we'll allow different reduced 
matrix elements for each. A purely conventional redefi- 
nition of the reduced matrix element is (J'HT^H J) = 

(-l)V2fc+T [£ fc i(l + (-l) k SP) T M k i(l - (-l) k 6P)] 
for the transverse case and (J'HT^HJ) 

-l) fc +V(2A; + l)C fe |(l + 



L (o)l 
-l) k SP) for 



the 



-P—(- 

longitudinal case 10 . SP is the product of initial 
and final meson parities and E(M) indicates an elec- 
tric (magnetic) multipole. After some manipulation 
we obtain the multipole decomposition of the helicity 
amplitudes, 



M(X, = ±) = J2 V 2JTT + (-i)"^) T M k \{l - {-l) k 5P)\ (fc T ; J'X ± l\JX) 



2fcH 


- 1 


2J- 


- 1 


2fcH 


- 1 


2J- 




- 1 



M(X, = 0) = J2 VItTT^ 1 + (-l) k 8P)(kO; J'X\JX). 



For the particular case under study, J = 0, J' = 1; this transition has one transverse multipole - E\ and one 
longitudinal multipole - C\ . 



M{X-, =±) = E x (q 2 ) 
M{X 7 = 0) = -C^q 2 ) 

Contracting the complex conjugate of (|A1|) with a photon polarisation vector yields, 

M(X, = ±) = (m 2 y - m 2 s )B + ( 9 2 ) - 9 2 S_(g 2 ) 

W 7 = 0) = [(m£ - m| + 9 2 )i?-(? 2 ) - (3m 2 , + m| - q 2 )B + (q 2 )} 

Solving (|A3|) . (|A4ft for B±{q 2 ) allows us to write the current matrix element as a multipole expansion. 



(A3) 



(A4) 



(S(ps)\rmV(Pv,r)) = fTV) E^q 2 ) n(q 2 )e»(p v ,r) - e(p v ,r).p s (p^PvPS ~ m 2 v p») 



i Ci(q 2 ) . _, . [ , .„ 2 u 2 u 

H T^mvqpi', rj.ps pv-Ps{pv + vsy ~ m sPv ~ m vPs 



(A5) 



r 



Note that this is expressed entirely in terms of invari- 
ants and covariant quantities and hence can be used in 
any frame. One can check that the tensor coefficients 
of the form-factors are orthogonal, indicating the inde- 
pendence of E\,C\. The invariant quantity £l(q 2 ) = 
{pv-Psf-m 2 v m 2 s = \[{m v -ms) 2 -q 2 ][{m v +m s ) 2 -q 2 ], 
and takes the simple value m||g| 2 in the rest frame of the 
decaying scalar. 

The decomposition is further constrained by 



Minkowski time-reversal invariance of the matrix 
element; this is a symmetry of the system if the meson 
states are stable. Time-reversal is implemented as an 
anti- unitary operation (T). As such the statement of 
invariance for the matrix element of a hermitian opera- 
tor is (/3\A\a) = (a\T AT _1 |/3) = (0\T AT _1 |a:))* 
where |<5) = T\a). States of definite momen- 
tum, spin and z-componcnt of spin transform as 
T\p,J,r) = ((— l)' 7 ~' r | — p,J,—r), where ( is an 
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arbitrary unphysical phase that we can choose indepen- 
dently for each particle tvpe[5q|. For the scalar and 
vector we choose £ = +1 - with this the vector has a 
real decay constant, fv defined in the standard way: 
(0^(0)7^(0) \V(q,r)) = m v f v e^(q,r). Note that to 
get a real pseudoscalar decay constant one needs to 
choose (p = — 1 (see the next appendix). 



The vector current transforms as Tj M T 1 = P%j u ■ Us- 
ing the relation e fJ, *(p, r) = ( — l) r+1 P£e I/ (— p, — r) one can 

show that Ei(q 2 ) and 1( ^- > are real. 

V<? 2 

An equivalent procedure for the axial- vector transition 
(1 + «-> 1~) yields the decomposition 



(Atf A ,r A )\f(0)\V(Py,ry)) = ^-^^(p^ -p v )„x 



E\ {q 2 ) (pa + Pv ) P ( 2mA [e* (pa , r A ) -Pv } ^ (Pv , ry ) + 2m v [e (p v ,r v ).p A ]el (pa , r A ) 



+ M 2 (q )(pa + pv)p\2m A [e*(p A ,r A ).p v }e v (pv,r v ) - 2m v [e(p v , r v ).p A ]<C(pA, r A ) 
I -4n(q 2 )e:(p A ,r A )e p (p v ,r v ) 



+ (pa +Pv) p (m A - my + q 2 )[e*(p A , r A )-Pv] e v (p v ,r v ) + (m\ -m v - q 2 )[e(p v ,r v ).p A ] el(p" A ,r A ) ) 



(A6) 



APPENDIX B: MINKOWSKI AND EUCLIDEAN 
N-POINT FUNCTIONS 

We obtained all our Lorentz decompositions in 
Minkowski space, but perform lattice computation in Eu- 
clidean; in this appendix we outline a simple way to effect 
the mapping between the two and describe some discrete 
symmetry constraints. 



1. Two-point functions 

In Minkowski space the following fermion bilinears are 
hcrmitian: 

S(x) = ip(x)ip(x) 
P(x) = ip(x)i'y 5 ip(x) 
V»(x) = ip(x)^ip(x) 

A"(x) =4)(x) 1 tl l 5 4>(x) 
T^(x) = i)(x)o^^( x ). (Bl) 

The transformation of spin-0 and spin-1 fields under 
time-reversal can be written 

with arbitrary phases £. The choices in IjBlll correspond 
to 



so that the decay constants defined in Q are real: 

wfp = (p^p(omo)Yi 5 m\p(-Q))r 
= -iPtpy f* P <p 

=> fp = fp, 

m v f v e»(q 7 r) = (P^-l) 1 ^ W(0) 7 ^(0)|U(-<f, -r)»" 
= Cv(-l) 1 - r m v f v P^*(-q,-r) 
=>• fv = fv- 

In addition, the overlap factors Z must be real: 

Z P (p) = (0|m 7 5 d(0,0)|F(p)) 

= ((0|r- 1 Tm 7 5 d(0, 0)T- l T\P(p))) * 

= ((0|(-l)m 7 5 d(0,0)Cp|P(-p)))* = Z* P (-p), 

since for the truly local operator, Z is not actually a func- 
tion of p and where for a rotationally-invariant smeared 
operator Z = Z(\p\). 

The mapping between Minkowski and Euclidean space- 
times is effected by the transformations, 

t -> -it 
. o 



74 

ilk 
-75, 



(B2) 



Cs 



-i; 



-i; Cv 



-i; 



where the tilded quantities are Euclidean (e.g. the 
gamma matrices satisfy { 7m , 7i /} = 28^). 
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Consider as an example the vector two-point function 
in Minkowski space, 



r^(p, t) = J2 e^V^W, t) [</¥V>(o, o)] 



(3 



E 



(4 



(AH 



■e * 



p 1 ^ 



where we inserted a complete set of states and used the 
Minkowski space completeness expression for polarisa- 



tion vectors Er eAI (^ r ) e * I/ (^ r ) = (-.9^ + £ ^~J- In 
the lattice calculation we actually compute a Euclidean 
quantity: 

f « = -(E^M^ 5 ^' ^ 5 ^^'^' ^} 

= ^6^(^(^,^(0,0)) 

= (-i) a rC 2 )«(s-it) 



E 



-e v> 



p l p J 



where with this derivation it is clear that all energies 
and momentum should be interpreted as the usual real 
Minkowski variants. All that we had to do was apply the 
mapping (|B2|) . Note that we only ever consider polari- 
sation vectors in Minkowski space, where they are easily 
defined. 



2. Three-point functions 



In Minkowski space, inserting two complete sets of 
states gives 



,-iEft f e -i{Ei-Bf)t 

2E f (p f ) 2EM) 



(0|#7V(0, 0)1/(29), 77)) (f(p f , r f )\m@, 0)\i(p h r,)} ((01^(0,0)^^:,^)) 



The Minkowski space-time Lorentz decomposition of the 
matrix element was discussed in the previous appendix. 
We actually compute the Euclidean variant, 

f %@f, Ql tf,t) = Cfcrcir {3)fri (pf,q; -itf, -it) 



where e.g. I\ = CiT 1 . As an explicit example consider 
the T] c form- factor where Tf = = ij 5 and r = 7 M . 
Then using 1B2|) we have 



T^ kp {p f ,q- -itf, -it) = lmtf\ kP {pf,q; if, t) 



r& p ~< 0p (p f ,?,-it f ,-ii) 



RcT 



(3) 

P74P 



(Pf,qjtf,t). 



One should perform the mapping on each transition com- 
puted to ensure that one extracts the correct complex 
component with the right sign. 



APPENDIX C: SCALE SETTING ON 
ANISOTROPIC LATTICES 



We outline how to set the scale of dimensionful quan- 
tities when one has differing spatial and temporal lattice 
spacings. The lattice action, written in terms of dimen- 
sionful quantities is S = J2x* o^atipQip. If we scale the 

w 3 /2 

fermion fields by ip = a s ip then we get a lattice spac- 
ing independent action, ipQip, if Q = o-tQ where the 
anisotropy will appear in the spatial derivative operator. 

On the anisotropic lattice a bosonic field would have a 
mass term in the action proportional to a^atm 2 ^ 2 , 
so that if we scale the mass as fa = atm (as we must 
since it appears in Euclidean correlators like e~ mt ), then 

/ „3 

we have to scale the boson field as 

If wc define the creation/annihilation operators by the 
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continuum decomposition 11 
d 3 p 1 



4>(x) 



(2tt) 3 2E P 



(a\p)e ipx +a(p)e-' lp - a 



and use pi = rii (3-momenta must scale with a s as 
they appear in the dimensionless combination p.x) we 
have 



1 



at 



L 3 a 3 2E ff 



a\p)e tp - x + a{p)e 



so that the dimensionless creation/annihilation operator 
is a — ^p^a. Then since a single particle state is de- 
fined by \N(p)) = cr\0), we define a dimensionless single 
particle state by \N(p}) = a t |0) = ^~^\N(p)). It is 

easy to check that the resolution of the identity takes the 
following form 



1 



E 

N 



d 3 p 1 

(2tt) 3 2Ep 



\N(p))(N(p)\ 



N 



2Ep 



which has no dependence on either lattice spacing. One 
can then show by insertion of the complete set of states 
that the two-point function calculated in a lattice simu- 
lation is, as expected, 



f ( 2 > = 4e" 
2E 



IE 



The scaling of Z is easily found 



Ze~ lp - S = (0|#>(aO \N(p)) 



1 



(0\^(x)\N(p)), 



so that the dimensionless Z that appears in the lattice 
simulation is Z = \J a 3 atZ. 

The pscudoscalar decay constant, /p, has mass dimen- 
sion 1 and is defined by 

(0|<K0) 7 'W(0)|P(<Z)) = ifpcf- 

Taking the temporal component and scaling to dimen- 
sionless quantities we have 



1 



(0|^(0) 7 U 7 5 ^(0)|P(<?)) =ifE, 



9"t ' 



so that fp = y^fp- Thus to obtain the physical decay 
constant from the lattice value we calculate 

/ = r 8/ V7. 

The same formula applies to the vector decay constant. 

As an example of setting the scale of a dimensionful 
transition form-factor factor, consider the scalar to vector 
transition: 

(S(p s )\4>^\V(pv,r)) = £ x (QV(fr,r) + ... 

The LHS scales to the dimensionless version leaving only 
one factor of a t in the denominator and hence E\ = atE\. 
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